Controlling double ionization of atoms in intense bichromatic laser pulses 
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We consider the classical dynamics of a two-electron system subjected to an intense bichromatic 
linearly polarized laser pulse. By varying the parameters of the field, such as the phase lag and the 
relative amplitude between the two colors of the field, we observe several trends from the statistical 
analysis of a large ensemble of trajectories initially in the ground state energy of the helium atom: 
High sensitivity of the sequential double ionization component, low sensitivity of the intensities where 
nonsequential double ionization occurs while the corresponding yields can vary drastically. All these 
trends hold irrespective of which parameter is varied: the phase lag or the relative amplitude. We 
rationalize these observations by an analysis of the phase space structures which drive the dynamics 
of this system and determine the extent of double ionization. These trends turn out to be mainly 
regulated by the dynamics of the inner electron. 

PACS numbers: 32.80.Rm, 05.45.Ac 
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I. INTRODUCTION 

Multiple ionization of atoms takes place at intermedi- 
ate and high intensities when they are subjected to laser 
pulses [1| . The routes by which these multiple ionizations 
occur can be quite complicated. Two main routes have 
been identified for the double ionization of two-active- 
electron systems (like the helium atom subjected to lin- 
early polarized laser pulses): a sequential process and a 
nonsequential one [l[. In the sequential process the field 
ionizes one electron after the other in an uncorrelated 
way. In contrast, the electron-electron correlation plays 
an active role in the ionization of the second electron 
in the nonsequential process. More specifically, the field 
picks up an electron, and sends this electron to the re- 
maining ion core to dislodge the second electron, accord- 
ing the so-called recollision (or 3-step) scenario @, H| . 

It has been observed experimentally in a broad range of 
intensities that non-sequential double ionization (NSDI) 
[also referred to as correlated double ionization (CDI) [J]] 
yields can be many orders of magnitude larger than in the 
sequential double ionization (SDI) [also referred to as un- 
correlated double ionization (UDI)] counterpart. Vary- 
ing the intensity gives rise to a characteristic "knee" 
shape for the double ionization yield as a function of 
the laser intensity. This strong effect due to electron- 
electron correlation has been observed experimentally 0- 
[l3j and also theoretically using quantum 11 1 [13-12(1 . s emi- 
classical (2ll - [23j and classical simulations [4 Il8l. l2ll. HJ- 

m. 

In this article, we investigate the dynamics resulting 
from the interaction between helium atoms and a bichro- 
matic linearly polarized laser field using a classical model. 
Classical means are entirely sufficient to produce the 
strong electron correlation needed to explain double ion- 
ization phenomena [25|. Atomic phenomena in intense 



perturbative regime of laser-atom interactions |30H42j . 
Although the classical mechanisms of double ionization 
remain qualitatively the same as for a monochromatic 
laser field [29( , the double ionization yields show a strong 
sensitivity to changes in the parameters of the laser field, 
namely the phase lag between the two colors and their 
relative amplitudes. We explain this dependence using 
an analysis in terms of phase space structures. Our main 
result is that all these processes are regulated primarily 
by the dynamics of the inner electron which is the elec- 
tron that remains bounded after the other one has been 
ionized. Remarkably, this electron, rather than being a 
spectator, is in fact the driver in the ionization process. 

We consider a classical model with one spatial dimen- 
sion for each electron which corresponds to the direction 
of polarization of the laser field. We denote the posi- 
tions of the two electrons by x and y and set the origin 
at the nucleus. We call the canonically conjugate mo- 
menta of x and y, p x and p v resp ectively. We consider 
soft Coulomb potentials 0, SB-ll [H, |H defined by two 
softening parameters a and b for the interaction between 
charged particles. The Hamiltonian of the system in the 
dipole approximation is obtained by the sum of the ki- 
netic energy of the electrons, the soft Coulomb potentials 
between the two electrons and between the electrons and 
the nucleus, and the interaction of the electrons with the 
field: 
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bichromatic fields can be considered as an extension of 
the principles of coherent control into the nonlinear non- 



\J (x — y) 2 + b 2 

— -, - + {x + y)E{u>t), (1) 

\Jx z + a z ^y z + a z 

where E (uit) is a bichromatic laser field, defined as: 

E (cot) = E Q f(ut) (sin (ut) + a sin (kut + 5)), (2) 

where Eq is the amplitude of the dominant mode, uj its 
frequency, a < 1 the ratio of the amplitudes of the two 
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colors of the laser field, k characterizes the mode lock- 
ing between the two components, and 8 is the phase 
lag. Since the second laser is an additional energy in- 
put to the system, the intensity is scaled with a as 
lo oc Eq (l + a 2 ). The system has two and a half de- 
grees of freedom, one for each electron and one half for 
time dependence. The envelope function / is chosen as 
a trapezoid with 2 laser cycle (I.e.) ramp up, 4 laser cy- 
cle plateau, 2 laser cycle ramp down, and zero elsewhere. 
In what follows, the wavelength is chosen to be 780 nm 
(corresponding to a frequency of w = 0.0584 a.u.) and 
k = 3. Except when we investigate the role of a, we set 
a = 0.45 for numerical simulations, which corresponds 
to an additional intensity input of 20%. We consider the 
initial state to be the ground state of He, i.e. 8 g = —2.24 
a.u. [25j, |44| . We choose a = 1 to prevent auto- ionization 
(so that the energy surface is bounded in phase space). 
The second softening parameter b is chosen so as to allow 
a significant energy exchange during a recollision. Here 
6 = 1@. 

For insight into the dynamics of typical trajectories 
and to interpret the mechanisms at play which explain 
the statistical analysis, we present two kinds of studies 
as functions of the parameters of the field, Eq, a and 8 in 
what follows: In Sec. HH a statistical analysis of a large 
ensemble of trajectories, and in Sec. lIIII an analysis of the 
phase space structures. 
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FIG. 1: Probability of double ionization as a function of the 
intensity for a = 0.45. The red crosses correspond to 8 = 0, 
the black circles to 8 — it/2, the green squares to 8 = n. We 
also indicate the range of variation of 1^ with 8 as given by 
Eq. (|14[) (vertical dashed-dotted lines). The vertical magenta 
continuous line labels the intensity chosen in Fig. [2j [Tl and [8] 



A. Varying the phase lag 8 



II. STATISTICAL ANALYSIS 



Starting from the classical model (JTj), we perform a sta- 
tistical analysis of a large ensemble of trajectories (typ- 
ically 150,000) initially on the ground state (with a mi- 
crocanonical distribution [29|]). At the end of the pulse, 
a trajectory can exhibit zero, one or two ionized elec- 
trons. In our calculations, an electron is considered ion- 
ized when its energy (defined as the sum of its kinetic 
energy and the Coulomb interaction with the nucleus) 
is positive, which corresponds to an electron free from 
the Coulomb potential well. A convenient observable is 
the probability of double ionization where the two elec- 
trons have a positive energy at the end of the pulse. This 
probability is commonly plotted as a function of the laser 
intensity. A signature of NSDI is a characteristic knee 
shape: In the intermediate intensity range, these yields 
are much higher (by several orders of magnitude) than 
expected from the uncorrelated sequential double ioniza- 
tion mechanism At high intensity, the yields reach 
saturation, corresponding to double ionization of all the 
electrons of the ensemble. In what follows, we investigate 
the double ionization yields versus the intensity when the 
phase lag 6 and the relative amplitude a are varied. 



In Fig. [T] we display several double ionization probabil- 
ity curves as a function of the laser intensity for different 
values of the phase lag 5 (for a — 0.45). 

In the intermediate range of intensities (up to satura- 
tion of double ionization), we observe a significant vari- 
ability of the double ionization probabilities and conse- 
quently of the characteristic features of the knee. More 
precisely, in the intermediate intensity range the level of 
double ionization yields greatly vary with S by a factor of 
3; however, when varying 5 the range of intensity where 
the knee is located (around 5 x 10 14 W • cm -2 , vertical 
dashed-dotted lines in Fig. [TJ does not change signifi- 
cantly. In the low intensity regime i.e., for intensities 
lower than 8 x 10 13 W • cm" 2 there is a very weak de- 
pendence of 8 in these yields. The same conclusions hold 
for other values of a. We note a significant impact of 
the phase lag in the resulting double ionization probabil- 
ity, meaning that adjusting this parameter can act as a 
control knob for both sequential and non-sequential pro- 
cesses. One advantage of varying the phase lag between 
the two components of the laser field is that no additional 
input of energy is required (as opposed to varying other 
parameters like a) since the mean value of the electro- 
magnetic energy does not depend on 6. 

In Fig. [2] we represent the double ionization yield as a 
function of S for the intensity labeled by the vertical ma- 
genta continuous line in Fig. [TJ The resulting curve has a 
bell shape with a maximum slightly shifted before 8 = n. 
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FIG. 2: Double ionization probability as a function of S for 
a — 0.45 and the intensity labeled by the vertical magenta 
continuous line in Fig. [T] (I = 4.78 x 10 14 W • cm" 2 ). 

More specifically, we observe a fast variation of the dou- 
ble ionization yield with changes of the phase lag around 
5 = and a flattening around 5 = tt. This confirms the 
high sensitivity of the double ionization probability to 
the parameter 5. 

Except for very low intensity (i.e. below 8 x 10 13 W • 
cm" 2 ) we observe qualitatively the same impact of S on 
the double ionization yield (a bell shape almost centered 
around 5 = tt). The precise value of the phase lag for 
which the probability is maximum depends on the laser 
parameters. 

The value of 5 for which we observe the first saturation 
of double ionization occurs when S = tt. This value of 
S does not depend on a (and is defined as soon as a is 
non-zero). We define Jw (5) the intensity, as a function 
of the phase lag, for which saturation of double ioniza- 
tion is reached. Then, at high intensity, we note a great 
variability of (5) (the differences are of the order of 
10 16 W-cm" 2 ) which is more than one order of magnitude 
larger that the range of variation of the knee. 

B. Varying the relative amplitude a 

After varying 5 at fixed a, we keep S fixed at tt in or- 
der to investigate en's influence on the double ionization 
probability. We chose 5 = tt because it corresponds to 
the phase lag where the peak field amplitude is maxi- 
mum (see Sec. IIIIDj) . In Fig. Owe represent the double 
ionization probability as a function of the intensity for 
different values of a. As for the variation of 5, we no- 
tice no significant dependence on the parameters at low 
intensities (below 8 x 10 13 W • cm" 2 ). In the intermedi- 



FIG. 3: Double ionization probability as a function of in- 
tensity for different values of a. Blue circles correspond 
to a = (single laser), Red crosses to a — 0.10, green squares 
to a = 0.45, black plus signs to a = 0.75 and magenta trian- 
gles correspond to a — 1.00. Here S = n for all curves. 



ate intensity range, we notice that the probability greatly 
varies with a (more than a factor two) whereas the in- 
tensity at which the characteristic knee is developed does 
not change significantly. In the meantime, the UDI com- 
ponent also increases with a, and we notice a decrease 
of the intensity 1^ where double ionization is saturated. 
As a consequence we notice a flattening in the knee, the 
two CDI and UDI curves merging. As a is increased, 
the second color of the laser field becomes more impor- 
tant, hence 5 is expected to have a more prominent role 
in controlling the system because the variation in proba- 
bility becomes larger. This is confirmed by Fig. U where 
we have represented double ionization probabilities as a 
function of the intensities for a = 0.1 and a = 0.75. 

In summary, the statistical analysis shows these three 
main trends when changing the parameters of the field: 

1. High sensitivity of the actual probability of double 
ionization, 

2. High sensitivity of the intensity of saturation of the 
double ionization yields, 

3. Low sensitivity of the intensity of knee regime (i.e. 
the CDI maximum). 

In the next section, we analyze the dynamics of Hamil- 
tonian (|TJ) in terms of phase space structures and we find 
results that support and explain these trends. 
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FIG. 4: Double ionization probability as a function of the 
intensity for different values of S. We keep a constant and 
equal to a = 0.1 and a = 0.75, respectively. 



FIG. 5: Typical nonsequential double ionization. We dis- 
play the position (upper panel) and energy (as defined in the 
text, lower panel) of each electron as functions of time. The 
parameters are I = 4.78 x 10 14 W-cm -2 , a = 0.45 and 5 = 0. 



III. NONLINEAR DYNAMICAL ANALYSIS 

For the dynamical analysis of Hamiltonian (JXJ) we fol- 
low the one performed in Ref. [2!| for the monochromatic 
laser pulse. The analysis proceeds in the following way: 
First, without laser field (Eq — 0), an outer electron and 
an inner one were identified from a periodic orbit analysis 
of the chaotic motion of typical trajectories. When the 
field is turned on, if the inner electron remains bounded, 
the atom may undergo a recollision which follows from 
the 3-step picture @, y] (the field picks up the outer elec- 
tron, and upon reversal of the field, brings it back to the 
core and the two electrons share the collision energy). 



A. Effective models 

The dynamics of both electrons are accurately repro- 
duced by reduced models except during recollisions [29j . 
These models allow one to derive the main properties of 
the double ionization curves versus intensity. We begin 
by investigating the behavior of typical trajectories. In 
Fig. [5j we represent the position of each electron (upper- 
panel) and their energy (kinetic energy plus soft Coulomb 
potential, lower panel) as functions of time for a typ- 
ical trajectory. We observe that one electron is first 
quickly ionized (the "outer" electron) while the other 
remains close to the core (the "inner" electron). Once 
ionized, the outer electron is driven mainly by the laser 
field and experiences large excursions from the core. In 
contrast, the inner electron experiences a competition be- 
tween Coulomb interaction with the nucleus and the laser 



excitation. Without loss of generality, we assign x to the 
outer electron and y to the inner one. In the sample 
trajectory of Fig. [5] the inner electron (red curves) does 
not ionize until it collides with the outer electron (blue 
curves). Double ionization occurs at approximately 2.5 
laser cycles, which can be confirmed by examining the 
energies of both electrons in the lower panel. We see 
that both electrons have positive energy after this time 
(and thus are ionized). From this qualitative analysis, 
two reduced effective Hamiltonians are identified, H out 
for the outer and Hi n for the inner one: 

2 

H out (x,p x ,t) = ^ + xE(cot), (3) 
H in (y, Py ,t) = gg- -^= +y E(ut), (4) 

where we have neglected the interaction with other 
charges for the outer electron and have allowed the in- 
ner electrons to feel both the laser field and the attrac- 
tion from the nucleus. For both electrons, we ignore 
the electron-electron interaction term since they are as- 
sumed to be far away from each other between two recol- 
lisions. The main advantage of the above effective models 
is that their dimensions are lower than the one of Hamil- 
tonian (they both have one and a half degrees of free- 
dom). Hamiltonian -ff ou t is even integrable in the sense 
that one can explicitly give the expression of the trajec- 
tories as functions of the initial conditions. Hamiltonian 
H m is not integrable but its dynamics can be analyzed 
qualitatively. In what follows, we analyze these two mod- 
els as the parameters of the field are varied and correlate 
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their dynamics with the double ionization curves plotted 
in Sec. [TT] in order to support and explain the observed 
trends. 



B. dynamics of the outer electron 

The trajectory of the outer electron [as governed by 
Eq. ([3J] is given by 

x(t) = xq + p t + — | (smut + — sin (kuit + 8)j (5) 

Eq /a \ 

Px (t) = Po H I cos cot + — cos (kujt + 6) I , (6) 

uj V k ) 

where xq + E§oik~ 2 u~ 2 sin (5 and po + 
i?o^ _1 (l + ak^ 1 cosS) are, respectively, the posi- 
tion and momentum of the outer electron at time t = 0. 
An important quantity is the maximum return kinetic 
energy of the outer electron: It gives the maximum 
amount of energy the outer electron can bring to the 
inner electron to trigger CDI. In addition, if the outer 
electron recombines, it gives the cutoff in harmonic 
generation for photon emission (46| . In order to compute 
this quantity we need to consider a shift of the field since 
the electron leaves the core at a certain phase 4> of the 
field. The effective Hamiltonian is transformed into 

H aat = ^+xE(ut + 4>). 

Maximizing the kinetic energy S max — p x (t r ) 2 /2 with re- 
spect to the recollision time t r and the initial phase <p 
leads to the familiar condition that the maximum return 
energy is associated with trajectories leaving the nucleus 
with zero momentum (46[. As a function of the param- 
eters of the laser field, this maximum kinetic energy is 
equal to 

£max = K (a, S) Up, (7) 

where U p — Eq/(Aoj 2 ) is the ponderomotive energy and 
k is only a function of a and 6. For a — 0, we have 
k(0,5) w 3.17 2]. In Fig. [SJ we represent n(a, S) as a 
function of 6 for various values of a. We notice that the 
amplitude of the variations of k with 6 increases with a. 
For a fixed value of a, the maximum return energy k (in 
units of Up) of the outer electron can vary significantly 
with 6 (up to more than 2). For all values of a, the 
maximum of k occurs at 6 ~ 1.04, while its minimum 
is obtained for S « 3.99. We notice that none of these 
values are related to the critical values of the phase lag for 
which we observe maximum of double ionization (around 
5 = 7r) and saturation of double ionization. Thus, the 
origin of the amount of double ionization does not solely 
rely on the dynamics of the outer electron. It has to be 
complemented by the dynamics of the inner electron. 




FIG. 6: Maximum return energy of the outer electron as 
a function of S in units of U p . Red squares correspond to 
a — 0.10, blue triangles to a = 0.45, and green circles to 
a — 1.00, while the dashed black line corresponds to the one 
laser system (a = 0). 



C. Dynamics of the inner electron 



After the ionization of the outer electron, the inner 
electron experiences the combined action of the field and 
of the attraction from the nucleus only. Its Hamiltonian 
is given by Eq. (j4]). The competition between these two 
interactions leads to two distinct behaviors related to two 
distinct regions. Figure [7] represents these two regions in 
phase space: More specifically, we plot the distance of 
the inner electron from the nucleus after 10 laser cycles 
in the space of initial conditions (yo,p y fi) at t = 0. Note 
that the color map is in logarithmic scale. The red and 
green curves are two different Poincare sections of an in- 
ner electron trajectory, taken stroboscopically with the 
laser period. They indicate that the dark region is filled 
with invariant tori where electrons are bound to the nu- 
cleus (for all time) . The light region corresponds to ioniz- 
ing electrons (after a transient time which depends on the 
initial conditions, see the darker structures which spiral 
out). The rationale behind this dynamics is as follows: 
If the electron is close to the nucleus, then the Coulomb 
interaction dominates over the field and the electron is 
bound. If it is sufficiently far, the field dominates and 
the dynamics of the electron is unbounded, in a similar 
way as the outer electron. The maximum distance the 
inner electron can be from the nucleus while remaining 
bound is a function of the laser parameters. 
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FIG. 7: Laminar plots of inner electron dynamics as given by Hamiltonian Q as a function of the initial conditions (yo,Py,a, t = 
0). For each initial condition, we represent the distance of the electron from the nucleus after 10 laser cycles. The distance is 
plotted in a logarithmic scale. In the left panel 5 = and in the right panel 5 — n. The dashed curves are Poincare sections 
of two typical trajectories in the bound region taken stroboscopically with the laser period. The parameters are a = 0.45 and 
I = 4.78 x 10 14 W ■ cm -2 in both plots. 



D. Double ionization probability 

Varying certain laser parameters leads to strong vari- 
ations in the size of the bound region, which is pivotal in 
explaining the probability of double ionization. Hamilto- 
nian (01 is not integrable. However, a close inspection of 
the dynamics of the inner electron in the bound region 
reveals that it is fairly insensitive to the field so we con- 
sider another effective Hamiltonian in the bound region, 
namely 



(8) 



which is integrable, its phase space being foliated by pe- 
riodic orbits. There exists, at least locally, a canonical 
transformation which maps Hamiltonian ((5J) into action 
and angle variables. The action is defined by 



A = — 
2 



When the field is present, most of the outer invariant tori 
are broken, whereas most of the inner ones close to the 
nucleus persist. We denote by A m (Eg, a, 8) the action 
of the last invariant rotational torus which roughly cor- 
responds to the boundary of the inner region. Figure [5] 
illustrates the dependence of A m (Eo,a,6) with 6 as a 
and Eq are kept fixed for an intensity in the UDI regime 
(/ = 4.78 x 10 15 W • cm -2 ). Numerically, we determine 
A m by computing the volume of phase space for which 



the inner electron never ionizes during 2 laser cycle ramp 
up and 4 laser cycle plateau (see Fig. [5]). At least qual- 
itatively, the same trend is observed in the whole UDI 
regime (as well as the CDI regime), namely an inverted 
bell shape with fast variations of A m around 6 = with 
a flattening and minimum centered around 5 = tt. 

An approximation to the size y m of the inner region can 
be obtained by finding the local maximum of the effective 
potential (soft Coulomb potential plus laser excitation) 
at peak field amplitude. It is implicitly given by 



M (a,6)E -- 
where M(a,5) is given by: 



(vl + 1) 



3/2 ' 



(9) 



M(a,S)= max Isinr + a sin (for + 6)\ , (10) 

rG[0,27r] 

and corresponds to the peak field amplitude of the laser. 
In Fig.[5]we display the peak field amplitude as a function 
of 5. In the bound region, the energy of the inner elec- 
tron, as defined by Eq. (j4]), does not vary much. Then, 
it is possible to relate it to the action A through the 
formula [4| 



H in (A) 



V2- 



where 7 — — 9\/2/16 (~ —0.8). It can be inverted to 
express the action explicitly as a function of the energy 



A =^ ln ( 1 + ^(^ + 2 ) 



(11) 
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FIG. 8: Action A m of the outermost invariant torus as a 
function of 5 for a = 0.45 and I = 4.78 x 10 14 W ■ cm -2 . The 
solid line corresponds to the statistical data, while the dashed 
dotted line corresponds to the theoretical model, as given by 
Eq. (HU). 
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FIG. 9: M (a = 0.45, 8) as a function of 8, as given by 
Eq. (HOJ). 



H m into Eq. (II]) 

where y m is solution of Eq. |H1 We compare this predic- 
tion with numerical evaluation of the action in Fig. [8l A 
careful analysis of A m , as given by Eq. (fl2"|) shows that 
A m is symmetric and minimum around S — tt for all a. 

Sequential double ionization (UDI) corresponds to the 
situation when one electron is first ionized by the field, 
leaving the inner electron in the unbound region, so that 
the field sequentially ionizes the remaining electron. The 
amount of UDI is related to the volume of the bound 
region, therefore to A m : the smaller the inner region 
(and thus A m ) the more likely it is for the inner electron 
to be found in the unbound region as well. Nonsequen- 
tial (CDI) double ionization corresponds to the situation 
when the inner electron is first left in the inner region, 
and is then propelled, through a recollision, in the un- 
bound region under the impact of the returning outer 
electron. As a result, it is easier for the outer electron 
to ionize the inner one if the bound region is small. Ei- 
ther way, for both CDI and UDI, we are expecting higher 
double ionization probability when A m is minimum (i.e. 
around 5 — n) and a strong sensitivity of the double ion- 
ization yield where A m varies fast (i.e. around 5 = 0). 
It corresponds to what is observed in Fig. [5] The small 
discrepancy between S — n and the actual maximum is 
due to nonlinear energy exchange processes during the 
recollision. In Fig. |51 we see that the kinetic energy of 
the outer electron upon return is not symmetric about 
5 = n. In particular we see that the maximum return en- 
ergy is rather small around n in comparison with slightly 
lower values of 6 which explains that the maximum is 
expected to be slightly smaller than tt. 

The rate at which A m changes with 6 provides an ex- 
planation for the sensitivity of the probability to changes 
in 6 which was discussed in Sec. Ill Al The action changes 
quickly in S in the regions around 5 = and slowly 
around 6 — tt, which is the same trend observed for the 
rate of change in probability with 5 from Fig. [5] This 
explains Observation [I] from the statistical analysis, in 
Sec. HH namely that the double ionization yield is highly 
sensitive to the choice of phase lag. 



E. Saturation of the double ionization probability 



The outermost position y m at peak field amplitude is 
associated with zero momentum, such that its energy H m 
is equal to 

2 

H m = H in (y m , 0, t m ) = - E M (a, 6) y m - 

VV™ + 1 

Finally, the critical action A m is obtained by plugging 



As explained before, CDI corresponds to the case 
where both electrons are left in the unbound region. An 
estimate of the size of the bound region is implicitly ob- 
tained through Eq. @. An approximation of the ampli- 
tude of the field at which there is no inner region leads to 
an estimate of the intensity at which double ionization is 
saturated, meaning that all the trajectories double ion- 
ize. Equation (0) has a real solution if EqM is smaller 
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than 4/(3 V3) which leads to 



3V3M(a,S)' 



The field is 
M (a, ?r) = 1 



maximum for 8 = 7r and its value is 
a. For this value of 8 the critical value of 
the intensity at which all the trajectories double ionize 
is minimum. This is what is observed in Fig. [T] and [4] 
A numerical estimate of the range of variation of this 
critical intensity as a function of 8, corresponding to the 
amplitude of variation of intensity between Eq (a, 0) and 
E^\a,Tr) gives 1.1 x 10 16 W • cm" 2 at a = 0.45 which 
is what is observed in Fig. [TJ A more precise approxi- 
mation of the value of the intensity (as a function of the 
laser parameters) at which the bound region disappears 
is obtained by following the elliptic periodic orbit (with 
the period of the laser field) at the center of the elliptic 
island [29]. It explains Observation [5] from the statistical 
analysis, in Sec. HH namely that the intensity of satu- 
rated double ionization greatly varies with the choice of 
the laser parameters. 




FIG. 10: Expansion of , where k and M have been solved 
for numerically. For this plot a — 0.45, k — 3, and cj = 
0.0584. 



F. Maximum of CDI 

When one electron is bound to the nucleus after the 
other one has ionized, the only way through which it 
leaves the nucleus is by an exchange of energy with the 
outer electron when this one comes back to the core 
upon sign reversal of the laser field. Through this recol- 
lision the inner electron jumps from one invariant torus 
to another one (depending on the amount of energy ex- 
changed), and it can even jump from one invariant torus 
to the unbound region. After it enters into the unbound 
region, it ionizes (solely due to the action of the field) . A 
nonsequential double ionization corresponds to the situ- 
ation when the exchange of energy allows the inner elec- 
tron to jump from the bound region to the unbound one, 
while the outer electron keeps enough energy to remain 
ionized. The amount of energy exchanged as well as the 
size of the bound region are the two main reasons for 
the sensitivity of the nonsequential double ionization to 
changes in the parameters a and S. 

Using an equal sharing of the energy of the inner 
and outer electron (29|, we derive a prediction for the 
intensity when the CDI probability is maximum. At 
{UiPy) — (0,0), when the inner electron is bounded, its 
energy is U- m — —2, while at (y,p y ) = (y m ,0), when the 
inner electron is on the boundary (on the invariant torus 
with action A m ), its energy is U ln = — 2/ y/ y 2 ~ + 1 . This 
leads to a difference in energy of A£ v = 2-2/ v / y 2 m + l 
which is the minimum energy required to ionize the inner 
electron in one recollision. This energy comes from the 
outer electron, and when 

AS V = (13) 



the outer electron gives enough energy to ionize the inner 
electron (in one recollision) while keeping itself ionized as 
well. 

The equal sharing relation can be solved in a manner 
very similar to the case of the one laser system [29[ by 
making the assumption that double ionization occurs at 
peak field amplitude, as given by Eq. ([9]). Combining 
Eq. (fTBl with Eq. © yields the following approximation 



^vH^ +£ #)l (14) 

For numerically computed k and M we can see from 
Fig. [TO] that 1^ does not vary significantly with 8. Over 
the entire range of 8, the critical intensity j( c ' varies by 
approximately 2 x 10 14 W-cm~ 2 which is compatible with 
the variations of 7^ as observed in Fig.[TJ This variation 
is small compared to the variations of the laser intensity 
from 10 13 to 10 17 W • cm" 2 over the entire CDI and UDI 
regimes. It is also small compared to the variations of 
/W which are of the order of 10 16 W • cm 2 . It confirms 
Observation [3] from the statistical analysis, in Sec. HH 
namely that the knee regime (and particularly the inten- 
sity for which CDI is maximum) is fairly insensitive to 
the laser parameters. 



IV. CONCLUSION 

Using the statistical and nonlinear dynamical analy- 
sis we showed that the laser parameters a and 8 can be 
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used to control both the probabilities of double ioniza- 
tion and the rates of change of these probabilities both in 
the CDI and UDI regimes. This control acts on the outer 
electron by regulating the amount of energy it can bring 
to the core during a recollision. For the inner electron, 
the control depends on its stability, by changing the size 
of the bound area. Those two modulations act together 
to determine the double ionization probability, but it is 
the inner electron which has the strongest impact on its 
variations. In a nutshell, using the inner and outer elec- 
trons reduced models, we were able to explain the three 
observations obtained with the statistical analysis of the 
system: Namely the high sensitivity of the double ion- 
ization yield and saturation of double ionization in the 



choice of the laser parameters and the low sensitivity of 
the knee regime in the same parameters. This control 
helps in regulating the proportion of sequential versus 
nonsequential processes which occur in this atomic sys- 
tem. 
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